Thermal equilibrium and efficient evaporation of an ultracold atom-molecule mixture 
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We derive the equilibrium conditions for a thermal atom-molecule mixture near a Feshbach res- 
onance. Under the assumption of low collisional loss, thermodynamical properties are calculated 
and compared to the measurements of a recent experiment on fermionic lithium experiment Q]. 
We discuss and evaluate possible collision mechanisms which can lead to atom-molecule conver- 
sion. Finally, we propose a novel evaporative cooling scheme to efficiently cool the molecules toward 
Bose-Einstein condensation. 
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I. INTRODUCTION 

Recent experiments witness a fast and remarkable 
progress in creating molecules from atomic Bose-Einstein 
condensates (BEC) @ 0, 0, degenerate Fermi gases 
[E IE or ultracold thermal gases 0, 0- These re- 
sults have initiated new pursuits toward molecular BEC, 
Cooper-paired Fermi gas [E [ljj as well as matter-wave 
interferometry based on ultracold molecules. 

In all these experiments, the creation of molecules 
is based on magnetically tunable Feshbach resonances 
[Til IT^ | , which allow interacting atom pairs to couple to 
molecules in a single internal quantum state. Extensive 
studies on atomic Feshbach resonances not only lead to a 
thorough understanding of the cold atom collision prop- 
erties, but also suggest intriguing possibility to convert 
atomic BEC or degenerate Fer mi ga s to a molecular BEC 
with high efficiency [3 El EE El . 

Converting atoms into molecules via a Feshbach res- 
onance is acc omp anied by other effects: three-body 
collisions 0, Il8l Il9j . atom- molecule interactions |y, 
molecule- molecule interactions pj , as well as many-body 
effects in quantum degenerate gases @- These processes 
can lead to a low conversion efficiency or a short molec- 
ular lifetime, both limit the ability to observe and ma- 
nipulate the molecules. As the detailed mechanisms of 
the conversion process require future investigations, it 
remains an experimental task to identify the best system 
and strategy to reach molecular condensation. 

In most recent experiments, molecules are created by 
ramping the magnetic field through a Feshbach reso- 
nance. High conversion efficiencies of 50% to 80% are 
found in degenerate Fer mi g ases [E IE 0> which agree 
with the calculations [2(J, |21| • Surprisingly, a recent ex- 
periment showed that a similar efficiency can be achieved 
in a thermal gas of fermionic b Li atoms at fixed magnetic 
field 0. Furthermore, the molecule sample can be pu- 
rified, trapped and remains stable with a lifetime of up 
to 10s. The results of Ref. [l| prompt us to investigate 
the thermodynamics of a non-degenerate atom-molecule 
mixture and the possibility to cool the molecules to a 
molecular BEC. 

In this paper, we introduce a thermodynamical model 
to calculate the atom-molecule conversion efficiency in 



a thermal gas (Sec. II). Using the model, we derive the 
atom-molecule (molecule-atom) conversion efficiency and 
compare the results with the experimental data (Sec. Ill 
and IV). To further cool the molecules, we propose a 
novel evaporative cooling scheme in an atom-molecule 
mixture (Sec. V). Finally, we identify the relevant con- 
version mechanisms and calculate the time scale of the 
proposed cooling process (Sec. VI). 



II. MODEL 

Based on the fermionic lithium system, we consider a 
classical thermal sample with N and N' atoms in two 
internal states and M molecules. The molecules indi- 
vidually consist of two non-identical atoms and are in a 
single internal state with constant binding energy E. We 
assume the system has negligible collision loss and energy 
exchange with the environment. 

The conservation of particle number and energy im- 
poses two constraints in the thermalization process: the 
former is simply N + M = const, and N' + M = const., 
while the latter is given by Ne n + N'e n > + M(e m — E) = 
const., where e n (e n <) and e m are the mean external en- 
ergy per atom and molecule. We assume the gas is dilute 
with negligible interaction energy among particles. 

The thermal equilibrium condition is obtained assum- 
ing that the atom-molecule mixture is in a quantum- 
mechanical canonical ensemble. This assumption is 
gen erally valid in the limit of large particle number 
|22| . Given the temperature T and single-particle par- 
tition function Z n (Z n /) and Z m c~ E l kT for atoms and 
molecules, we derive the equilibrium condition by min- 
imizing the free energy F — —kT hiZ, subjected to 
only the particle number conservation constraint, where 
Z '= Z^Z% 'Z^e- ME / kT (N\N'\M\)- 1 is the partition 
function of the system, and k Boltzmann's constant. The 
result constitutes the key equation we investigate in this 
paper, 



= <t>N4>N'e 



E/kT 



(1) 



Here 4>m = N/Z n (<j>M = M/Z m ) is the final atomic 
(molecular) phase-space density in the motional ground 
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FIG. 1: Thermal equilibrium in an atom-molecule mixture. 
Molecular states \v m ) are introduced to the two-body states 
of non-identical atom pairs \u,u'), where v, v' (f m ) are the 
motional quantum numbers of a single atom (molecule) . Pop- 
ulation in |0, 0) is the product of that in single-particle state 
\v — 0) and \v' = 0), namely, (f>N<j>N' ■ Molecular phase- 
space density is enhanced by the Boltzmann factor e E ^ kT 
when E > 0. 



state. 

This equation is valid for systems with atoms in two 
internal states and molecules in one state, as well as 
those with two atomic species, say, Rb atoms, Cs atoms 
with RbCs molecules. The only assumptions are the low 
collision loss and interaction energy. For single compo- 
nent systems with TV atoms and M molecules, Eq. Q is 
rewritten as </>m = 4>N eE ^ kT H- 

The thcrmalization condition given in Eq. can be 
understood in a simple picture. As the atomic phase- 
space density 4>n and 4>n' are the ground state popula- 
tions of the two atomic components, 4>n4>N' is the popu- 
lation of the non-identical atom pairs in the lowest two- 
atom state. Eq. (JJJ states that the motional ground state 
of the molecule is regarded as one of the two-atom states. 
An additional Boltzmann factor e E l kT accounts for the 
binding energy. This is illustrated in Fig. ^ Notice that 
this result is independent of the details and complexity 
of the molecule formation and dissociation mechanisms. 

Interesting consequences can be directly observed from 
Eq. l|T)l. At constant temperature, molecular density is 
proportional to atomic density squared, even in the case 
when the molecules are formed, say, by a three-body pro- 
cess, whose rate depends on density cubed. This is be- 
cause the reverse molecule dissociation rate also depends 
on atomic density. The balance of both the formation 
and dissociation processes in thermal equilibrium gives 
the density squared dependence. We will come back to 
the dissociation process in Sec. VI. 

In the following sections, we calculate the atom 
(molecule) number N (M) and temperature T in ther- 
mal equilibrium based on Eq. JQl and the conservation 
laws: we first determine the constants of the conservation 
laws from the initial conditions, then express partition 
function Z n , Z n i and Z m in terms of trapping poten- 



tial and temperature T; mean energy in terms of parti- 
tion function and temperature, = —Z^ x dpZi, where 
(3 = (kT)^ 1 . In cases where no analytic solutions are 
available, we solve the equations numerically. 

A special case is calculated here when the molecules 
and two atomic components are identically trapped Z n = 
Z n > = Z m with initial atom numbers Nq, Nq and 
molecule number Mq — 0. The conversion efficiency 
/ = 2M/(N Q + iV£) can be derived from Eq. Q and 
particle conservation law as 

/=(l + ^^e-Wf. (2) 
V Z<Pn<Pn> J 

When the atomic gas is tuned right on resonance E = 
with final phase-space density <pN = <Pn' = 1, an conver- 
sion efficiency of 50% is obtained. Notably, the unity 
phase-space density assumption is on the boarder of the 
applicability of our model. 

On the other hand, given an atom sample with low 
initial phase-space densities 4>n — 4>n^ *C 1 and zero 
binding energy, the conversion efficiency is simply / = 
4*n — 4*N' ~ 4>n - This result indicates that creating 
molecules from a cloud of thermal atoms by tuning the 
field right on Feshbach resonance is inefficient. In the 
following sections, we show that an appreciable gain in 
conversion fraction can be obtained when we tune the 
molecular state below continuum. 



III. CONVERTING ATOMS TO MOLECULES 

The atom-molecule conversion can be qualitatively un- 
derstood from Eq. QJ. When the molecular state is far 
below (far above) the atomic continuum, all particles 
should accumulate in the lower molecular (atomic) state 
and the atom (molecule) number is exponentially sup- 
pressed. This naive picture, however, is incorrect when 
one tries to convert cold atoms into molecules. At large 
binding energy, the internal energy released during the 
molecule formation process heats up the sample signif- 
icantly and reduces the atomic phase-space density. In 
the following calculation, we show that the final temper- 
ature goes up approximately linearly with the binding 
energy and the molecule fraction is therefore limited. In 
the limit of infinite binding energy, counter-intuitively, 
no molecules are formed. 

To show the conversion in the vicinity of a Feshbach 
resonance, we consider both atoms and molecules are har- 
monically trapped with identical single-particle partition 
function, Z n = Z n > = Z m = ^(1 - e - fia "/ fcT )- 1 , where 
u>i is the trap vibration frequency in the ith direction. 
The external energy of the particles is e„ = e' n = e m — 
3kT. The assumption that both species have the same 
trap vibration frequencies is generally valid for atoms and 
long-range atomic dimers in a deep far-detuned dipole 
trap where the trap depth and mass for a molecule are 
both twice as large as those for an atom. 
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binding energy E/k [uK] 

FIG. 2: Atom and molecule fraction in equilibrium starting 
with a pure atom sample. Trap parameters are specified in 
the text. In the upper figure, the atom fraction N/No (dotted 
line) and the molecule fraction M/No (solid line) are calcu- 
lated and compared to the measurements (solid circles) . The 
atom(molecule) fraction of 50% at large binding energy is ac- 
cidental. In the lower figure, the final phase-space density for 
molecules (j>M (solid line) and for atoms <f>M (dotted line) are 
shown together with the temperature T (dashed line). 




10"* 10' 1(T 10 10" 

initial atomic phase-space density <j> N 

FIG. 3: Maximum conversion efficiencies in molecular phase- 
space density 4>m and in molecule number M from a pure 
atom sample. The calculation is performed assuming an 
isotropic harmonic trap, an initial atomic phase-space den- 
sity <pN and atom number No = Nq. Notice that maximum 
molecular phase-space densities <j>M (solid line) and numbers 
M (dashed line) are obtained at different binding energies E. 



factor of 5 ~ 20 lower than (j)^ , while the maximum 
conversion fraction varies from 30% ~ 55%, shown in 
Fig. [3] This very weak dependence on the initial atomic 
phase-space density is remarkable and suggests that an 
efficient conversion of atoms into molecules in thermal 
clouds is possible. 



Based on the experiment parameters given in Ref. [l|: 
axial (radial) trap frequency of uj z (ui r ) = 2ir x 260 (0.39) 
kHz, initial atom number Nq = = 1.25 xlO 6 , molecule 
number Mo = 0, and temperature To = 2.5/aK, we cal- 
culate the molecule (atom) fraction, phase-space density 
and temperature in thermal equilibrium, shown in Fig. [21 
A maximum conversion efficiency of 2M/(Nq + Nq) = 
M/Nq ~ 50% agrees with the experimental results |24j. 

Remarkably, the maximum conversion efficiency is 
reached at the binding energy of E > 20/dK. This value 
greatly exceeds the initial temperature of the sample and 
qualitatively explains the displacement of the maximum 
recombination loss position by ~ 150G relative to the 
Feshbach resonance (2^,|2^|. However, a generalization 
of our treatment to quantum-degenerate gas is necessary 
to predict the exact held shifts in those experiments. For 
even more negative binding energy, an abrupt turn-off of 
the molecule formation rate is observed and no experi- 
mental data on the molecule number in thermal equilib- 
rium is available. 

To calculate the best atom-molecule conversion effi- 
ciency, we assume initially Ao = Nq atoms are in an 
isotropic harmonic trap with Z n — Z n i — Z m , phase- 
space density 4>n — <f>N' an d binding energy E, we cal- 
culate the maximum molecular phase-space density 4>m 
or conversion fraction M/Nq by varying the binding en- 
ergy. Over the range of 4>n = 10~ 4 ~ 1, the maximum 
phase-space density of the molecules <pM is typically a 



IV. CONVERTING MOLECULES TO ATOMS 

Converting a pure molecule sample into atoms, as 
demonstrated in Ref. |l|, shows different features. First, 
at large negative binding energy, all molecules dissoci- 
ate into atoms. Beginning with Mq = 3 x 10 5 molecules 
at 2.5/zK with the same trap parameters described ear- 
lier, we calculate the molecule (atom) fraction M/Mq 
(N/Mq) in thermal equilibrium, shown in the upper fig- 
ure of Fig. A full conversion from molecules to atoms is 
achieved when the molecular state is high above the scat- 
tering continuum. Shown tog ether with the calculation 
is the experimental data |2J|. 

The agreement between the experiment and calcula- 
tion is excellent as there are no free parameters in the 
calculation. Both the experiment and calculation sup- 
port the possibility to dissociate molecules into atoms 
even when the binding energy is positive and large com- 
pared to the initial molecular temperature. This is due 
to the much larger phase space of the atomic scattering 
continuum as compared to that of the molecular bound 
state. 

Furthermore, a gain in phase-space density after ther- 
malization is predicted at small positive binding energies 
where temperature drops and phase-space density peaks 
up, shown in the lower hgure of Fig. 01 This cooling 
occurs due to two different processes. First, the dissoci- 
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FIG. 4: Atomic and molecular fractions in equilibrium start- 
ing with a pure sample of Mo molecules. In the upper figure, 
the fractions of molecules (solid line) and atoms (dotted line) 
are calculated and compared to the experimental data: atoms 
(solid circles) and molecules (open circles). In the lower fig- 
ure, the phase-space density of the molecules (solid line) and 
atoms (dotted line) are shown together with the final temper- 
ature (dashed line) 



ation process is endoergic and reduces the total external 
energy. Second, the total particle number, N + N' + M, 
increases after molecules dissociate and further reduces 
the mean energy per particle. When the equilibrium is 
reached, we find a gain in molecular (atomic) phase-space 
density of 2.4 (4.2), compared to that of the initial molec- 
ular phase-space density. This cooling effect is studied in 
more detail in Section V. 



V. NOVEL EVAPORATIVE COOLING IN AN 
ATOM-MOLECULE MIXTURE 

In conventional evaporative cooling schemes, energetic 
particles are removed from the trap. The remaining par- 
ticles rethermalize and acquire a lower temperature and 
higher phase-space density. 

The molecule-atom conversion process resembles the 
above process. By locating the atomic scattering con- 
tinuum above the molecular state with positive binding 
energy, molecules with higher thermal energy are more 
probable to collisionally dissociate into atoms. While 
these particles are immediately removed in the conven- 
tional evaporative cooling recipe, the atoms from the dis- 
sociated molecules have a greatly reduced thermal en- 
ergy and a rethermalization among these atoms and the 
remaining molecules is advantageous to cool the sam- 
ple. To evaluate the cooling performance, we consider 
an initially pure sample of Mo molecules in an isotropic 



harmonic trap with phase-space density <f>M i tempera- 
ture To and binding energy E — r/kTo, where r\ defines a 
truncation parameter, in analogy of the energy cutoff in 
the conventional evaporative cooling scheme. After the 
equilibrium is reached, we calculate the mean evaporative 
cooling efficiency based on 7 = — ln((j)M / 4»m ) /\ n (M / M ) 
for a wide range of truncation parameter 77, shown in 
Fig. where 4>m (M) is the molecular phase-space den- 
sity (number) in equilibrium. At constant 77, we find the 
atom-molecule thermalization permits a better cooling 
efficiency 7 than does the conventional one. 

On this basis, we propose a scheme to continuously 
cool molecules toward high phase-space density and pro- 
vide a quantitative estimate on its performance. Given 
a harmonically trapped molecule sample with negligible 
collisions loss, temperature T and total external energy 
K, we control the molecular binding energy according 
to the cloud temperature and truncation parameter ac- 
cording to E = rjkT. We assume atoms resulting from 
molecule dissociation are quickly thermalized with the 
sample before they are removed. The removal of only 
atoms from an atom-molecule mixture is demonstrated 
in Ref. 0. During the thermalization process, a small 
molecular fraction e that dissociates will lower the total 
external energy by 5K = —Keri/3; next, the thermaliza- 
tion process shares the energy among the remaining 1 — e 
molecules and 2e atoms and lowers the temperature by 
ST = —Te(f] + 3)/3. Given the phase-space density in 
a harmonic trap <f>M = MZ^ 1 oc MT~ 3 , we obtain the 
evaporation efficiency as 7 = rj+2. The 77+2 factor comes 
from the increase in particle number, together with the 
endoergic nature of the dissociation process. Compared 
to the conventional evaporation with -f c = 77 + k — 4 |27| , 
where < k < 1, the new scheme with, say, 77 = 4 can 
have an evaporative cooling performance comparable to 
that in conventional method with 10 > 77 > 9, shown in 
Fig. 

However, this evaporation method involves additional 
atom-molecule conversion and thermalization steps and 
can potentially be slow due to either a slow thermal- 
ization rate or atom-molecule conversion rate. In the 
two-component fermion system, in particular, a large s- 
wave scattering length and elastic collision cross section 
for atom-molecule (molecule- molecule) scattering near an 
atom-atom Feshbach resonance are expected 0, . A 
fast thermalization rate on the order of (lms)" 1 is pos- 
sible from recent experiments. The conversion rate from 
molecules to atoms, however, depends on the details of 
the atom-molecule interaction mechanism and will be dis- 
cussed in Sec. VI. 



VI. REACTION MECHANISMS AND 
EVAPORATION RATE 

To estimate the conversion rate or evaporation rate, 
specific molecule formation and dissociation mechanism 
need to be identified. We will consider Li system as a 
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FIG. 5: Evaporation efficiency in an atom-molecule mixture. 
Evaporation efficiency 7 during the thermalization process is 
evaluated for r\ = E/kTo = 10 -4 (dashed line), and 77 = 10° 
(dotted line). Estimations based on the proposal (solid line) 
and on conventional evaporation (dash-dotted line) are also 
shown for comparison. No collision loss is considered in the 
calculation. 

specific example. 

In a dilute gas with negative binding energy, molecules 
dissociate according to the following process, 

Li 2 <-> Li + Li'. (3) 

This spontaneous process is dominant since it is den- 
sity independent. The molecular dissociation rate in this 
case is characterized by the Feshbach resonance linewidth 
|29j |. In typical experiments, typical Feshbach resonance 
linewidth of 0.1MHz ~ 10MHz indicates a very short dis- 
sociation time of < 1/j.s. 

When the binding energy is tuned to positive values 
E > 0, the process described by Eq. © is forbidden due 
to the energy and momentum conservation. The lead- 
ing binary processes to dissociate molecules and the cor- 
responding reversed ternary or quaternary processes to 
form molecules are 

Li 2 + Li *-> Li + Li' + Li, (4) 
Li 2 + Li 2 <-> Li 2 + Li + Li' (5) 
Li 2 +Li 2 «-> Li + Li' + Li + Li', (6) 

Notice that the dissociation (formation) processes are 
shown from left to right (right to left) and are allowed 
only when the total energy of the incident channel is suffi- 
cient to support the internal energy of the outgoing chan- 
nel. 

In thermal equilibrium, each molecule formation rate 
is exactly balanced by the corresponding dissociation 
rate, due to detailed b alancing . For the above pro- 
cesses, we have C\mn = i?in 2 n', C 2 m 2 = R^mnn' and 
C^m 2 = R 3 n 2 n' 2 , respectively, where Cj (R4) are the 
associated dissociation (formation) rate coefficients, n, 



n' (m) are the atomic (molecular) densities and X de- 
notes the averaged value of X in a canonical ensemble. 
To directly relate the dissociation and formation coef- 
ficients, we assume n (n') atoms and m molecules are 
uniformly distributed in a box with unity volume. The 
single particle partition function is given by Z n — Z' — 
2~ 3 / 2 Z m = A^J, where \ dB = h(2nm kTy 1 / 2 is the 
thermal de Broglie wavelength of the atom, tuq atomic 
mass and h Plank's constant. Using Eq. QJ, we obtain 
the relationship between the coefficients in processes de- 
scribed by Eq. (QJ, Eq. © and Eq. ®: 

Ci = R 1 h- 3 (nm kT) 3 ^e~ E / kT , (7) 
C 2 = R 2 h- 3 (irm kT) s / 2 e- E / kT , (8) 
C 3 = R 3 h- e (7rm kT) 3 e- 2E / kT . (9) 

The above "formation-dissociation" relation reveals the 
rate coefficient when the counterpart coefficient is calcu- 
lated or measured. 

In samples which consist of mostly atoms, the domi- 
nant atom- molecule conversion process is given in Eq. 
The formation reaction is widely studied in cold atoms 
systems and is called three-body recombination 0, 0, 
19]. For a two-component fermionic system, the rate co- 
efficient is i?i = 167a e kT/h for E/kT > 1 19]. Us- 
ing the relationship E — h 2 /moa 2 , and Eq. Q, we 
immediately obtain the dissociation coefficient C± = 
2>.7hh 2 m~ 3/2 {kTfl 2 E- 3 e- E / kT for E/kT > 1 HJ. 

The situation is different when we start with a pure 
sample of molecules, since the process in Eq. Q does 
not happen in the absence of atoms. Comparing Eq. (JSJ 
and Eq. JjJJl, we identify Eq. JSJ as the dominant disso- 
ciation process at low temperature since C3 is exponen- 
tially suppressed relative to C 2 . This suppression can be 
understood as the molecule-molecule collision energy in 
Eq. |JBJ| should be sufficiently high to support four atoms 
in the continuum, while in Eq. JSJ, only two atoms are 
in the continuum. 

In contrast to that in Eq. (QJ, the formation process 
in Eq. JSJ involves scattering of three non-identical par- 
ticles: Li 2 , Li and Li' and is expected to scale as i? 2 ~ e° 
at low collision energy e in the three-body scattering 
channel. Consequently, dissociation coefficient scales as 
C 2 ~ e 2 . Compared to R\ ~ e and C\ ~ e 3 [3(J , the dom- 
inant conversion process in an atom-molecule mixture at 
low temperature limit is actually Eq. Q, which involves 
four atoms. A detailed calculation will be necessary to 
quantitatively determine either C 2 or i? 2 . 

Knowing the possible dissociation mechanism for 
molecules, we estimate the speed of the proposed evap- 
orative cooling based on Eq. ©, in which the rate co- 
efficients are known. Given the evaporation parameter 
of 77 = E/kT = 4, temperature of 2.5/iK, and a small 
atom fraction with density n = 2 x 10 n cm~ 3 , we obtain 
an evaporation rate of G\n= (2.5s) -1 , which is indeed 
much slower than the two-body collision rate. However, 
in conjunction with the predicted evaporation efficiency 
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of 7 ~ 5, we expect an increase of molecular phase-space 
density by 5 orders of magnitude can be achieved in 6s, 
which is still within the lifetime of the molecular cloud 
10s reported in Ref. Q. Furthermore, during the evap- 
oration process, temperature decreases and the dissocia- 
tion process in Eq. Q will eventually dominate and speed 
up the evaporation. 

VII. CONCLUSION 

We have described the equilibrium condition and con- 
version mechanisms in a thermal mixture of atoms and 
molecules under the assumption of negligible collisional 
loss. Near a Feshbach resonance, the molecular state 
extends the phase space of that of two atoms. A full 
thermalization in this extended phase space provides the 
equilibrium condition given in Eq. Q and permits a 
quantitative estimation on the thermodynamical prop- 
erties of an atom-molecule mixture. 

In particular, we calculate the atom-molecule conver- 
sion efficiencies under various conditions and the results 
agree with the experimental data very well P, |24| . We 
show that by properly locating the binding energy of 
the molecules, the conversion of a significant fraction of 
atoms to molecules is possible even in a thermal gas with 
low phase-space density. Our result suggests an alterna- 
tive scheme to reach quantum degeneracy: first, convert 
the fermionic atoms into bosonic molecules; second, evap- 
oratively cool them to a molecular condensate and finally, 



convert them back to fermionic atoms. The last process 
would allow a creation of atomic Cooper pairs 0, El • 

To cool the molecules to a molecular condensate, we 
suggest a novel cooling method based on the endother- 
mic molecule-atom conversion process at positive binding 
energy. We estimate a much higher evaporative cooling 
efficiency can be obtained as compared to conventional 
cooling in magnetic traps or dipole traps. 

To estimate the time scale of thermalization and atom- 
molecule conversion, we identify a new collision process 
given in Eq. (JSJ , which dominates the atom- molecule con- 
version at low temperature. The evaporation rate of the 
proposed scheme is estimated and an increase of molecu- 
lar phase-space density by more than five orders of mag- 
nitude within the typical lifetime observed in the Li2 
thermal gas is expected pj. This result highlights great 
prospect of attaining a molecular Bosc-Einstein conden- 
sate. 
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